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In many models of CP violation, the electric dipole moment (EDM) of a heavy charged 
or neutral lepton could be very large. We present an explicit model in which a heavy 



O 
(N 
O 

ix. 

^ ■ neutrino EDM can be as large as 10 e-cm, or even a factor of ten larger if fine-tuning is 

allowed, and use an effective field theory argument to show that this result is fairly robust. 
■ We then look at the production cross section for these neutrinos, and by rederiving the 

Bethe-Block formula, show that they could leave an ionization track. It is then noted that 
the first signature of heavy neutrinos with a large EDM would come from e^e~ NN'y, 
leading to a very large rate for single photon plus missing energy events, and the rate 
and angular distribution are found. Finally, we look at some astrophysical consequences, 
including whether these neutrinos could constitute the UHE cosmic rays and whether 
their decays in the early universe could generate a net lepton asymmetry. 
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1 Introduction 



Most of the unknown parameters of the Standard Model come from fermion masses and 
mixing angles. They are put in by hand, and we have no real understanding of their 
origins. Many models that try to explain their values exist, most involving additional 
symmetries, but more experimental data will be needed to distinguish between them. 
Much of the excitement concerning neutrino masses and mixing angles is caused by the 
hope that their values will facilitate our understanding of the flavor problem. 

A potentially valuable source of information about the flavor problem may come from 
electric and magnetic dipole moments. Just as the Yukawa couplings form a matrix in 
generation space, the interaction of two fermions with a photon will also be a matrix 
in generation space. The real and imaginary parts of the diagonal elements will lead to 
the magnetic and electric dipole moments; the off-diagonal elements wil lead to flavor- 
changing radiative decays. If one can measure these moments, one would obtain valuable 
hints of physics beyond the Standard Model and would learn more about the origin of 
flavor. 

This paper is concerned with the electric dipole moments (EDMs) of leptons. There 
is a possibility of substantial improvement in the experimental bounds on EDMs in the 
near future. Proposals exist [|T| to lower the current bound[^ on the muon's EDM by six 
orders of magnitude, and substantial improvement in the bound on the electron EDM[^ is 
possible. Combining limits on the weak dipole moment of the tau with U{1) symmetry^ 
improves the current bound by a factor of thirty. 

How big might one expect the EDMs to be? In the Standard Model, they are negligi- 
bly smallQ. However, they can be much larger in extensions of the Standard Model. In 
multi-Higgs models, the EDM of the muon can easily be as large as 10~^^ e-cm, within 
reach of planned experiments^. In leptoquark models, the muon and tau EDM's are 
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typically 10~^^ e-cm and 10~^^ e-cm, respectively [||. Left-right models]^ have a muon 
EDM which is typically of the order of 10~^^ e-cm, and in the minimal supersymmetric 
standard model(MSSM) |]10[ , the electron EDM is above the experimental bounds if the 
phases are all unity. The point is that a wide variety of models give EDMs that can be 
observed in the near future. 



Babu, Barr and Dorsner||ll|| have discussed how the EDMs of leptons scale with the 
lepton masses. In many models, such as the Standard Model and the MSSM, they scale 
linearly with the mass. However, in a large number of models, such as multi-Higgs, 
leptoquark and some flavor symmetry models, the EDM scales as the cube of the lepton 
mass. In these models, the tau EDM will be 5000 times larger than that of the muon 
(it should be noted that the electron EDM in some of these models receives a two-loop 
contribution which varies only linearly, and thus it need not be negligible). 

Given cubic scaling, the EDM of a heavy lepton (charged or neutral) could be quite 
large. For example, a 100 GeV heavy lepton would have an EDM a billion times larger 
than the muon's. If the latter is in the expected range of 10~^^ e-cm, then such a heavy 
lepton would have an EDM as large as 0.01 e-fermi. This EDM would then dominate the 
electromagnetic interactions of these leptons, drastically changing their phenomenology. 
In particular, if a heavy neutrino acquired such an EDM, it could even leave an observable 
ionization track. 

Although originally motivated by cubic scaling models, the possibility of a heavy 
lepton having a huge EDM is worth studying in a model-independent way. We will see 
below that plausible models exist in which a heavy neutrino has a very large EDM, even 
though these models do not have cubic scaling. In this Article, we will concentrate on the 
possibility that a heavy neutrino could have a large EDM. In an earlier paper by one of 
us[0], the differential production cross section for heavy leptons (charged and neutral) 
with large EDMs was studied, and the fact that heavy neutrinos with large EDMs could 



leave an ionization track was discussed (although no explicit models were mentioned). 
This Article extends this earlier work substantially. In Section II, we explicitly present a 
model in which the EDM of a heavy neutrino can be of O(10~^^) e-cm. This will serve as 
an existence proof that relatively simple models can exist in which this occurs. Current 
bounds on such large EDMs will be noted. In Section III, the cross section for heavy 
neutrino production is presented, and the Bcthc-Block formula for a large-EDM heavy 
neutrino traversing matter will be calculated. This extends slightly previous work. The 
result will be that an ionization track is observable, but one would need to either modify 
existing detectors or construct a new detector. Since it is unlikely that this would occur 
without some preliminary evidence that such a large-EDM neutrino exists, we look, in 
Section IV, at the process of initial and final state radiation, e^e~ — > NN^, which would 
give an enormous rate for single-photon plus missing energy events, and would provide 
impetus for looking explicitly for large-EDM heavy neutrinos. Finally, in Section V, 
we discuss whether these heavy neutrinos could be candidates for the ultra-high energy 
cosmic rays, and also consider the implications for baryogenesis through leptogenesis, 
and in Section VI, we present our conclusions. 

2 Plausible Models 

Are there plausible models in which a huge EDM for a heavy neutrino occurs? If one 
does not rely on any particular model, then a simple effective field theory argument can 
be used. Suppose one assumes that CP violation is due to some sort of new physics at 
the TeV scale. Then one can write the effective low-energy, dimension-five Lagrangian 
as 

C = —LLd^yi'^^LRF^y. (1) 
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If A is 0(1) TeV, and the unknown coefficient is of 0(1), then this yields a very large 
EDM of approximately 10^^^ e-cm for the lepton. 

Of course, in realistic models, one expects the EDM to occur via a loop, suppressing 
the unknown coefficient, and there may not even be a dimension-five operator. What 
about specific models currently in the literature? A model with cubic scaling is that of 
Bernreuther, Schroder and Pham[|l^], in which CP violation occurs in the Higgs sector. 
In this model, the EDM of a heavy lepton (they consider the top quark, but the results 
are unchanged) is constrained by the electron EDM, and the maximum value is 0(10^^^) 
e-cm. In other models, such as the model of Babu, Barr and Dorsncr[|TT| in which their 



parameter c = 0, the electron EDM does not give any such constraints and a large EDM 
is allowed. 

Here, we will provide an explicit example of a model with a large EDM. Consider a 
model containing singlet fermions, -El,_r and Nl^r and two charged singlet scalar fields, 
hi. The fermions can couple to the charged scalars, /ij, z=l, 2, in the following form 

Y.^rNlK + Z.^lNrK + h.c. (2) 

where Yi and Zi are complex Yukawa couplings and can generally include flavor indices. 
The hi can also mix with each other and can be diagonalized by a complex unitary matrix 
U: hi = UijHj, with mass eigenstates Hi. Two such charged singlets are necessary so that 
a phase redefinition cannot eliminate the CP violation. For simplicity in presentation, 
we will look at the single H case, and just assume that the phase redefinition does not 
occur. In this case, the coupling is 

E{S - Pi^)NH + h.c. (3) 

where S and P are arbitrary complex couplings. 

It is straightforward to calculate the EDM due to the above interaction. The one-loop 
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Figure 1: One- loop feynman diagrams contributing to a large EDM for a heavy neutrino. 

feynman diagrams are listed in Fig. 0. The EDM is defined at = 0, and is given by 

2_n^_„^i^r^.oo*^ f\,. {l-z){l + 2z) 



dN{q' = 0) = e-^Im{SP*) / dz — '-j- -^^ ^ (4) 

which, numerically, in units of e-cm, is 

\ Me J 47r 7o -^(1 - 2)a2 + (1 - 2) + v ; 

where a = Mn/Me and 6 = Mh/Me- 

The integral is plotted in Fig. We see that EDM's of 10~^^ e-cm are quite possible, 
especially in the region of paramter-space in which the E is heavier than the N or H, 
if one chooses ^"^^^J" ^ ~ 1, and can even be an order of magnitude higher if one allows 
for some fine-tuning. We have not considered the unlikely region of parameter space in 
which Mn > Me + Mh, since the width of the in this case is comparable to its mass, 
and the definition of the EDM is no longer straightforward. 

We thus have presented a model in which the EDM of a heavy neutrino is very 
large, of 0(10^^^) e-cm. While this model is, admittedly, designed to have such a large 
EDM, its ingredients are very simple, and one can imagine that it may be part of a 
more complicated model. The main point is that having a heavy neutrino with a large 
EDM is certainly not excluded, and in this paper we are exploring the phenomenological 
implications. 
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Figure 2: Value of the integral in the expression for the EDM. 

Are there any current phenomenological bounds? If the heavy neutrino and charged 
lepton are part of a fourth chiral fermion family (albeit with a right-handed neutrino), 
then the S and T parameters will cause severe constraints, but if they are part of a 
vectorlike multiplet, there will be no such constraints. Very explicit models may be 
constrained, but there are no general bounds. One possible concern would be g-2 of 
the muon. Here, there will be an effect of a large EDM on the photon propagator. At 
the lowest order, the feynman diagram is illustrated in Fig. ^. As in the calculation of 
hadronic contributions to {g — 2)^, we calcualte the effect of a large EDM on {g — 2)^ by 
the dispersion integral [l^ 



^r""' = / , dsKis)ais),+,-^NN, (6) 

where K(s) is given in Ref. and a{s)e+e--.NN is given in Ref. [|T^], as well as 

below. We cutoff the upper bound of the integral. Letting m^r = 100 GeV, a^^^ = 
5.9 X 10"" (1.3 X 10"^°) for a cutoff of 1 (10) TeV. This is much smaller than the 



7 



Figure 3: Large EDM contribute to 
theoretical and experimental uncertainties, and is thus negligible. 

3 Direct detection 

The most dramatic effect of a large EDM of a heavy neutrino will be in the production 
cross-section and angular distribution. Escribano and MassoQ noted that the relevant 
f/(l) invariant effective operator is given by L^a^'^i'-y^LjiB^^, where B^^, is the U{1) field 
tensor. This gives a coupling to the photon, which we define to be the EDM, as well as a 
coupling to the Z which is the EDM times ta.n 6w- We will include this coupling to the 
Z. It turns out that the contribution from the Z has very little effect on the numerical 
results. In general, one should also include an operator coupling to the SU (2) field tensor, 
leading to a very different value for the CP- violating coupling to the Z. Rather than deal 
with two parameters, however, we just assume that the latter operator is smaller. If this 
assumption is false, then it will just make the cross-section even bigger, unless there is 
fine-tuning. Also, unless the CP-violating coupling to the Z is surprisingly large, it will 
have very little effect on the results. 

A discussion of the differential cross section for a heavy charged lepton can be found 



in Ref. II^. Here we are interested in heavy neutrino production. The differential 



8 



0.2 




Figure 4: Differential cross section for heavy neutrino production for various EDMs, in 
units of e-cm, for a lepton mass of 100 GeV. 



cross-section is given by 



da a" "^MU 1 , (l-4sin^g^)tangH^ „ \ 

dQAs\l s \ 8 sm* 26^;^ sm'' 29w J 

where 

2 .2^^^ 



Ai = D'ssin'e(l 



s 



4M^ 4M^ 

^2 = 1 + cos^^ sm^e + 8Cvcose + D^stan'^9w(sm'^0 + (l + cos^^)) 

s s 

A3 = 4D2s(sin2^ + ^!l^(l + cos2^)) (8) 

s 

where we have dropped the numerically negligible Cy terms, for simplicity. 

The differential and total cross-sections are given in Figures 4 and 5, for a heavy 
neutrino mass of 100 GeV. We see that for D = 0, the usual l+cos'^6+C cos 6 distribution 
for a lepton (where the cos 6' term is due to 7 — Z interference) is found. For a very 
large EDM, the distribution is completely dominated by the sin^ 6 contribution from the 
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Figure 5: The total cross section for heavy neutrino production for various EDMs, in 
units of e-cm, for a lepton mass of 100 GeV. 



electric dipole term. The angular distribution begins to deviate significantly from the 
D = distribution for D greater than 3 x 10^^^, and by -D = 3 x 10^^^ the distribution 
is very close to sin^ 6. One can see that EDM's of the size noted earlier will dramatically 



alter the angular distribution ||I2 



Note that since the cross-section varies as D^, for an EDM as ridiculously large as 1.0 
e-fermi, it would be almost a microbarn!! Cross sections this large will violate unitarity. 
The unitarity limit can be estimated by setting the relevant effective interaction strength, 
aDy/s, equal to unity. For y/s = 500 GeV, D = 10~^^ e-cm, and a ^ 1/125, this effective 
interaction strength is 0.25. Thus the unitarity limit will not be reached unless the EDM 
is larger than 10~^^ e-cm, and so the behavior discussed above will appear for a wide range 
of parameter space without violating the unitarity bound. Another way of saying this is 
to note that the larger the EDM, the smaller the scale at which the physics responsible 
for the effective interaction sets in, and for an EDM larger than 10~^^ e-cm, that scale is 
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less than y/s. 

Production of a fermion with a large EDM has been considered elsewhere. It is given 



in the context of a top quark EDM by Bernreuther et al.[13|, and given in the context of 



tau-pair production in Ref. ||T5|. This latter paper noted how one can use CP-odd angular 



correlations to search for a tau EDM, and this method has been used by experimentalists. 
However, there has not been any discussion of the possibility of a large EDM for heavy 
neutrinos, and here we see a unique and unusual signature. 

What is the purpose of this calculation? After all, it would not seem remotely possible 
to detect heavy neutrinos directly. This calculation would only be meaningful if the heavy 
neutrino were heavier than the charged lepton, and could thus be detected via its decay. 
Here, one could again look for CP-odd correlations, as discussed in Ref. [|l^]. However, 
most models have the charged lepton heavier than the neutrino, and thus the decay can 
only occur through mixing with the very light neutrinos. As discussed in detail in Ref. 
[ p^ , this mixing could be very small, and these heavy neutrinos could be effectively 
stable. As we will see below, however, it may be possible to detect these neutrinos 
initially through final state radiation, and then directly. 

Because the heavy neutrino has such a large electric dipole moment, it might be 
possible to detect it directly through its ionization loss. To estimate whether or not 
this might be possible, we follow the derivation of the classical Bohr formula outlined 
in Jackson||17||, replacing the electric field from a charge with the electric field from a 
dipole. This is a classical calculation; if one is not too far above threshold, the neutrino 
will be moving non-relativistically. Consider a heavy neutrino moving in the x-direction, 
and an electron at an impact parameter y = b. The impulse given to the electron is 
= I^cxi ^Edt. The result will depend on the orientation of the dipole. Suppose 
that the dipole is in the z-direction, transverse to the plane of the particle motion and 
the electron. Then the only nonzero electric field component is in the 2;-direction, and 
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Ez — -^-{h"^ -\-vH'^)~^/'^ , where the time of closest approach is defined to be i = 0, v is the 
velocity, and eD is the size of the electric dipole moment. Integrating, one finds that the 
impulse is |^^- Suppose that the dipole is in the y-direction. Then the electric field 
components are — ^^^^.s (3 sin 6> cos 6*) and Ey — ^^^^.s (3 cos^ 0—2), where — b'^+vH'^ 
and tan^ = b/vt. Integrating, the impulse in the x-direction vanishes, as expected by 
symmetry, and the impulse in the |/-direction is also |^;J2- Finally, if the dipole is 
in the x-direction, the electric field components are Ey — ^^^^3 (3 sin 9 cos 9) and E^ — 
(3 cos^ 9 — 1). Both of these integrate to zero, so there is no net momentum transfer 
in this case. Since the orientation is generally arbitrary, and we are only interested in 
a rough order-of- magnitude estimate, we will take the impulse to be |Ap| = |^^- 
Note that the impulse from an electric charge is just ^^f^^^ , so this result is expected 
dimensionally. This impulse is then converted into an energy transfer, AE = |Ap|^/2m. 
Jackson notes that the maximum energy transfer is AE^ax — and thus the 

minimum impact parameter is = D / {rwyv^) . The energy loss is obtained by 
cylindrically integrating over the impact parameter 



dE 
dx 



27rNZ / AE{b)bdb (9) 



and we obtain 



which is the corresponding formula to the Bohr result. Since we are doing a classical 
calculation, one can set 7 = 1. Note that the logarithm in the usual Bethe-Block for- 
mula is absent and the electron mass and particle velocity drop out. This is due to the 
extra power of b in the expression for the impulse, which arises on dimensional grounds. 
Plugging in numbers, one finds that ^ = lO^^D MeVg~^cm^. 

For EDM's of 10~^^ e-cm, this gives an energy loss which is approximately 10~^ of 
the usual energy loss for a charged particle. This is challenging for experimenters, but 
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not impossible, requiring specialized detectors. For example, if the heavy neutrino is 
produced in the decay of a charged lepton, then it would be produced in coincidence 
with a real or virtual W, which would eliminate backgrounds. Note that the event rate, 
for D = 10~^^ cm, is huge, corresponding to 100, 000 events for an integrated luminosity 
of 100 fb^^. This is several events per hour, and thus a specialized detector located on the 
outside of the main detectors could suffice, without needing a full 47r coverage (note that 
the iV's will emerge back-to-back). More importantly, however, is that initial and final 
state radiation will be huge, and thus standard searches for single-photon plus missing 
energy events will discover a huge signal, which could then be followed by a specialized 
search for direct detection. We now turn to the possibility of detection through single 
photon events. 

4 Detection through single photon events 

In the standard model, the process e^e^ vV^ has been calculated and used to de- 



termine the number of light neutrinos [|18], |T9[. The signature of the process is a single 



photon from initial state radiation plus missing energy. The process didn't rule out the 
existence of heavy neutrinos. Similarly, the process e^e~ NN'j can be used to study 
the effects of large EDM of heavy neutrinos. Now the single photon can also come from 
final state radiation because of the large EDM. 

The interaction of heavy neutrinos with the photon is —ieN Da j^ip'-y^q^ N , where q is 
the momentum transfer. The EDM is defined to be eD{q^ = 0). In this case, the photon 
momenta at the relevant vertices are different, and are certainly not zero. However, 
we presume that D{q^) does not vary too rapidly with g^, as is the case in which the 
momentum transfer is not substantially greater than the masses of the particles in the 
loop. The feynman diagrams are listed in Fig. The contribution due to the exchange 
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Figure 6: Feynman diagrams which contribute to the process e~^e NN'j. 



of the Z boson is neglected here since it is numerically small 1 20 1. 

The calculation is straightforward. It is clear that the rate is the sum of terms 
proportional to D^, and D'^. The contributions proportional to , and D'^ are 
presented separately in Fig. |^ for D = 10~^^ e-cm. For different D's, these curves scale 
appropriately. One can see that for D < 6 x 10^^^ e-cm, expected in realistic models, 
the term dominates. Because of this scaling, instead of falling as the cross section 
reaches a constant (it would grow as s if the D"^ term dominated. As discussed above, 
however, it does not reach the level at which unitarity is a big concern. The total cross 
section is then given, to a good approximation, by the term in Fig. 0. 

The cross section will be singular as the | cos 6*^1 tends to 1. We have used cuts in the 
photon energy and in the angle: ^ > 0.1 and | cos^^^| < 0.94. For example, at y/s = 500 
GeV, the cross section is about 1.5 x 10"^ nb (and scales as D^). This is a very large 
event rate, corresponding to about an event every few minutes (for D = 10~^^ e-cm) at 
a linear collider with a luminosity of 3 x lO^^cm^^sec"^. 

The angular distribution is also important. It is shown in Fig. ^ for D = 10~^^ e-cm, 

14 




10-^ ^ 

200 250 300 350 400 450 500 

(GeV) 

Figure 7: Different contributions to the total cross section for the process e^e~ —>■ NNj 
with the heavy neutrino mass of 100 GeV and D = 10~^^ e cm. For a positive EDM, the 
term is negative; the absolute value is shown. 
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Figure 8: Angular distribution of the cross section with the heavy neutrino mass of 100 
GeV and D = IQ-^^ e cm, and s^/^ = 590 GeV. 

and the resuhing forward-backward asymmetry is shown as a function of D in Fig. 
The angular distribution of the Standard Model is nearly forward-backward symmetric, 
but in this case there is a sizable asymmetry. For much of the parameter space, the 
asymmetry is —0.48. 

Clearly, the first signature of an EDM for a heavy neutrino will be a very large event 
rate for single photon events, with the forward-backward asymmetry measurement giving 
a confirmation of this particular interpretation of such a large rate. 

5 Astrophysical and Cosmological Significance 

One of the greatest mysteries in cosmic ray astronomy concerns the existence of cosmic 
rays with energies above the GZK bound[^. A substantial number of cosmic rays with 
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Figure 9: The forward-backward asymmetry as a function of the EDM. 
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energies well above 10^° eV have been observed, interacting in the upper atmosphere. 
Yet a charged particle's mean free path through cosmic microwave background, at these 
energies, is only about 50 megaparsecs, and no sources of such high energy particles 
are known to exist within this distance. Neutral particles will not have this problem. 
However, the fact that the cosmic rays shower high up in the atmosphere implies that 
the cross-section must be of strong interaction strength. The only strongly interacting 
long-lived neutral particle is the neutron, yet its lifetime is too short. Neutrinos will not 
interact high in the atmosphere, if they have conventional weak interactions. 

Several papers have considered the possibility that neutrino interactions might 
become strong at very high energies, due to exchange of Kaluza-Klein gravitons|22, 23|. 



One might wonder whether heavy neutrinos with a sufficiently large EDM might interact 
strongly enough in the upper atmosphere, but not interact strongly with the microwave 
background. Are these candidates for the ultra-high energy cosmic rays? One can redo 
the calcuation described above concerning ionization loss in a detector, but now looking 
at the ultrarelativistic limit. This is a tedious, but straightforward calculation. 

In this case, the cross section can not be high enough. Recall that the EDM is defined 
at = 0. However, here one would expect to be enormous. Certainly, in the first 
few interactions, it will be much larger than the mass-squared of the particles in the loop 
that produces a large EDM in the first place. As a result, the effective EDM will be much 
smaller, and the interaction will not be strong. These heavy neutrinos will not interact 
high in the atmosphere, and may not interact much more than conventional high energy 
neutrinos. 

Are there any other potential astrophysical or cosmological effects? One possibility 
concerns the generation of the baryon asymmetry. Suppose the heavy neutrino decays 
radiatively into light neutrinos. Then the huge EDM, which violates CP, may cause the 
decay rate into neutrinos to be different than the decay rate into antineutrinos. This 



would result in a large lepton asymmetry. Since non-perturbative splialeron interactions 
can convert a lepton asymmetry into a baryon asymmetry, this would result in a baryon 
asymmetry. The reader is referred to the article of Riotto and Trodden |^ for a nice 
review. This possibility is currently under investigation. 

6 Conclusions 

In this Article, we have explored the possibility that a heavy neutrino could have a huge 
electric dipole moment. It has been shown that plausible models exist that would lead 
to a dipole moment as large as 10~^^ e-cm, and that this would dramatically alter the 
electromagnetic properties of these neutrinos. The first signature of such a neutrino 
would come from a large enhancement of the single photon plus missing energy event 
rate. Then, one could actually look for the ionization track of these neutrinos in specially 
designed detectors. Cross sections and angular distributions are calculated as a function 
of the electric dipole moment. While these neutrinos can not explain the ultra high 
energy cosmic rays, they may have interesting implications for baryogenesis. 
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